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Although the rocket-based X-ray Quantum Calorimetry (XQC) experiment was designed for X- 
ray spectroscopy, the minimal shielding of its calorimeters, its low atmospheric overburden, and its 
low-threshold detectors make it among the most sensitive instruments for detecting or constraining 
strong interactions between dark matter particles and baryons. We use Monte Carlo simulations to 
obtain the precise limits the XQC experiment places on spin-independent interactions between dark 
matter and baryons, improving upon earlier analytical estimates. We find that the XQC experiment 
rules out a wide range of nucleon-scattering cross sections centered around one barn for dark matter 
particles with masses between 0.01 and 10 5 GeV. Our analysis also provides new constraints on 
cases where only a fraction of the dark matter strongly interacts with baryons. 

PACS numbers: 95.35.+d, 12.60.-i, 29.40.Vj 
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I. INTRODUCTION 

From Vera Rubin's discovery that the rotation curves 
of galaxies remain level to radii much greater than pre- 
dicted by Keplerian dynamics Q to the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) measurement of 
the cosmic microwave background (CMB) temperature 
anisotropy power spectrum [2| , observations indicate that 
the luminous matter we see is only a fraction of the 
mass in the Universe. The three-year WMAP CMB 
anisotropy spectrum is best-fit by a cosmological model 
with f2 m = 0.241 ± 0.034 and a baryon density that is 
less than one fifth of the total mass density. The cold 
collisionless dark matter (CCDM) model has emerged 
as the predominant paradigm for discussing the missing 
mass problem. The dark matter is assumed to consist 
of non-relativistic, non-baryonic, weakly interacting par- 
ticles, often referred to as Weakly Interacting Massive 
Particles (WIMPs). 

Although the CCDM model successfully predicts ob- 
served features of large-scale structure at scales greater 
than one megaparsec [3j , there are indications that it may 
fail to match observations on smaller scales. Numerical 
simulations of CCDM halos 0, i, i, 0, i, |, E3, El El 
imply that CCDM halos have a density profile that in- 
creases sharply at small radii {p ~ r~ 12 according to 
Ref. [ll). These predictions conflict with lensing ob- 
servations of clusters [H, E3 that indicate the presence 
of constant-density cores. X-ray observations of clusters 
have found cores in some clusters, although density cusps 
have also been observed [H, EE El- On smaller scales, 
observations of dwarf and low-surface-brightness galaxies 
PJ, El, S3, SI SI, SI Hi indicate that these dark mat- 
ter halos have constant-density cores with lower densities 
than predicted by numerical simulations. Observations 
also indicate that cores are predominant in spiral galaxies 



as well, including the Milky Way [25|, |26|, |27| . Numerical 
simulations of CCDM halos also predict more satellite 
halos than are observed in the Local Group [28|, |29( and 
fossil groups [30l |. 

Astrophysical explanations for the discord between the 
density profiles predicted by CCDM simulations and ob- 
servations have been proposed: for instance, dynamical 
friction may transform density cusps into cores in the 
inner regions of clusters [H , and the triaxiality of galac- 
tic halos may mask the true nature of their inner den- 
sity profiles |32j. There are also models of substructure 
formation that explain the observed paucity of satellite 

haios HH Si, m^. 

Another possible explanation for the apparent failure 
of the CCDM model to describe the observed features 
of dark matter halos is that dark matter particles scat- 
ter strongly off one another. The discrepancies between 
observations and the CCDM model are alleviated if one 
introduces a dark matter self-interaction that is compa- 
rable in strength to the interaction cross section between 
neutral baryons [H SI : 

odd 



= 8 x 10~ 25 - 1 x 10~ 23 cm 2 GeV-\ (1) 



where odd is the cross section for scattering between dark 
matter particles and TOdm is the mass of the dark mat- 
ter particle. Numerical simulations have shown that in- 
troducing dark matter self-interactions within this range 
reduces the central slope of the halo density profile and 
reduces the central densities of halo cores, in addition to 
destroying the extra substructure [13, S3] • 

The numerical coincidence between this dark mat- 
ter self-interaction cross section and the known strong- 
interaction cross section for neutron-neutron or neutron- 
proton scattering has reinvigorated interest in the pos- 
sibility that dark matter interacts with itself and with 
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baryons through the strong nuclear force. We refer to 
dark matter of this type as "strongly interacting dark 
matter" where "strong" refers specifically to the strong 
nuclear force. Strongly interacting dark matter candi- 
dates include the dibaryon [3l|, H2] , the Q-ball [43j] , and 
O-helium 0. 

Surprisingly, the possibility that the dark matter may 
be strongly interacting is not ruled out. While there are 
numerous experiments searching for WIMPs, they are 
largely insensitive to dark matter that interacts strongly 
with baryons. The reason is that WIMP searches are typ- 
ically conducted at or below ground level based on the 
fact that WIMPs can easily penetrate the atmosphere 
or the Earth, whereas strongly interacting dark matter is 
multiply scattered and thcrmalized by the time it reaches 
ground level and its thermal kinetic energy is too small 
to produce detectable collisions with baryons in WIMP 
detectors. Consequently, there are few experiments ca- 
pable of detecting strongly interacting dark matter di- 
rectly. Starkman et al. [45[ summarized the constraints 
on strongly interacting dark matter from experiments 
prior to 1990, and these constraints were later refined 
HI, 51) S3) S3 • The strength of dark matter interactions 
with baryons may also be constrained by galactic dynam- 
ics [45[, cosmic rays l45l l49j . Big Bang nucleosynthesis 
(BBN) jH, the CMB [5(|, and large-scale structure gtj. 

The X-ray Quantum Calorimetry (XQC) project 
launched a rocket-mounted micro-calorimeter array in 
1999 (IJ. At altitudes above 165 km, the XQC detector 
collected data for a little less than two minutes. Although 
its primary purpose was X-ray spectroscopy, the limited 
amount of shielding in front of the calorimeters and the 
low atmospheric overburden makes the XQC experiment 
a sensitive detector of strongly interacting dark matter. 

In this article, we present a new numerical analysis 
of the constraints on spin-independent interactions be- 
tween dark matter particles and baryons from the XQC 
experiment using Monte Carlo simulations of dark mat- 
ter particles interacting with the XQC detector and the 
atmosphere above it. Our work is a significant improve- 
ment upon the earlier analytic estimates presented by 
some of us in Refs. [H, HH because it accurately models 
the dark matter particle's interactions with the atmo- 
sphere and the XQC instrument. Our calculation here 
also supersedes the analytic estimate by Zaharijas and 
Farrar [52| because they only considered a small por- 
tion of the XQC data and did not include multiple scat- 
tering events nor the overburden of the XQC detector. 
We restrict our analysis to spin-independent interactions 
because the XQC calorimeters are not highly sensitive 
to spin-dependent interactions. Only a small fraction of 
the target nuclei in the calorimeters have non-zero spin; 
consequently, the bound on spin-dependent interactions 
between baryons and dark matter from the XQC experi- 
ment is about four orders of magnitude weaker than the 
bound on spin- independent interactions [52l |. 

This article is organized as follows. In Section |TT] we 
summarize the specifications of the XQC detector. We 



then review dark matter detection theory in Section [Ml 
This Section includes a discussion of coherent versus in- 
coherent scattering and how we account for the loss of 
coherence in our analysis. A complete description of our 
analysis follows in IIVI and our results are presented in 
Section |Vl Finally, in Section IVII we summarize our 
findings and compare the constraints to strongly inter- 
acting dark matter from the XQC experiment to those 
from other experiments. 



II. THE XQC EXPERIMENT 

Calorimetry is the use of temperature deviations to 
measure changes in the internal energy of a material. By 
drastically reducing the specific heat of the absorbing 
material, the use of cryogenics in calorimetry allows the 
absorbing object to have a macroscopic volume and still 
be sensitive to minute changes in energy. These detectors 
are sensitive enough to register the energy deposited by a 
single photon or particle and gave birth to the technique 
of "quantum calorimetry," the thermal measurement of 
energy quanta. 

The quantum calorimetry experiment [5l| we use to 
constrain interactions between dark matter particles and 
baryons is the second rocket-born experiment in the XQC 
(X-ray Quantum Calorimetry) Project, a joint under- 
taking of the University of Wisconsin and the Goddard 
Space Flight Center [HI HI- It launched on March 28, 
1999 and collected about 100 seconds of data at alti- 
tudes between 165 and 225 km above the Earth's surface. 
The detector consisted of thirty-four quantum calorime- 
ters operating at a temperature of 0.06 K; for detailed 
information on the XQC detector functions, please refer 
to Refs. [U These detectors were separated from 
the exterior of the rocket by five thin filter panes [5l| . 
The small atmospheric overburden at this altitude and 
the minimal amount of shielding in front of the calorime- 
ters makes this experiment a promising probe of strongly 
interacting dark matter. 

The absorbers in the XQC calorimeters are composed 
of a thin film of HgTe (0.96 fim thick) deposited on a sili- 
con (Si) substrate that is 14 /an thick. The absorbers rest 
on silicon spacers and silicon pixel bodies. Figs. [1] and 
[5] show side and top views of the detectors with the di- 
mensions of each layer. Temperature changes in all four 
components are measured by the calorimeter's internal 
thermometer. The calorimeters report the average tem- 
perature over an integration time of 7 ms in order to 
reduce the effect of random temperature fluctuations on 
the measurement. Multiple scatterings by a dark matter 
particle will register as a single event because the time it 
takes the dark matter particle to make its way through 
the calorimeter is small compared to the integration time. 

The detector array consists of two rows of detec- 
tors, with seventeen active calorimeters and one inactive 
calorimeter in each row, and is located at the bottom of 
a conical detector chamber. Within a 32-degree angle 
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FIG. 1: A vertical cross section of an XQC calorimeter. The 
relative thicknesses of the layers are drawn to scale, as are 
their relative lengths, but the two scales are not the same. To 
facilitate the display of the layers, the vertical dimension has 
been stretched relative to the horizontal dimension. 
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FIG. 2: A top view of an XQC calorimeter. The absorber is 
the top layer and underneath it lies the spacer, followed by 
the pixel body. These dimensions are drawn to scale. 

from the detector normal, the incoming particles only 
pass through the aforementioned filters. The five filters 
are located 2 mm, 6 mm, 9 mm, 11 mm and 28 mm above 
the detectors. Each filter consists of a thin layer of alu- 
minum (150 A) supported on a parylene (CH) substrate 
(1380 A). The pressure inside the chamber is less than 
10~ 6 Torr. At this level of evacuation, a dark matter 
particle with a mass of 10 6 GeV and a baryon interac- 
tion cross section of 10 6 barns, would have less than a 
20% chance of colliding with an air atom in the cham- 
ber. Therefore, we assume that the chamber is a perfect 
vacuum in our analysis. 

While the atmospheric pressure at the altitudes at 
which the detector operated is about 10~ 8 times the at- 
mospheric pressure at sea level, the atmospheric over- 
burden of the XQC detector is still sufficient to scat- 
ter incoming strongly interacting dark matter particles. 
Simulating a dark matter particle's path through the 
atmosphere requires number-density profiles for all the 
molecules in the atmosphere. These profiles were ob- 
tained using the MSIS-E-90 model 1 for the time (1999 



1 Available at http://modelweb.gsfc.nasa.gov/models/msis.html 
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FIG. 3: The points depict the MSIS-E-90 density profiles for 
the seven most prevalent constituents of the atmosphere above 
the XQC detector, and the lines show the piecewise exponen- 
tial fits used in our analysis. 



March 28 9:00 UT) and location (White Sands Missile 
Range, New Mexico) of the XQC rocket launch. 

During the data collection period, the average altitude 
of the XQC rocket was 201.747 km. At this altitude and 
above, the primary constituents of the atmosphere are 
molecular and atomic oxygen, molecular and atomic ni- 
trogen, helium, atomic hydrogen, and argon. The MSIS- 
E-90 model provides tables of the number densities of 
each of these seven chemical species. In our analysis, 
computational efficiency demanded that we fit analytic 
functions to these data. We found exponential fits for 
the density profiles in three altitude ranges: 200-300 km, 
300-500 km and 500-1000 km. The error in the proba- 
bility of a collision between a dark matter particle and 
an element of the atmosphere introduced by using these 
fits instead of the original data is 0.02%. Fig. [3] shows 
the number density profiles provided by the MSIS-E-90 
model and the exponential fits used to model the data. 

The XQC detector collected data for a total of 150 
seconds. During these 150 seconds of activity, the thirty- 
four individual calorimeters were not all operational at 
all times. Furthermore, events that could not be ac- 
curately measured by the calorimeters and events at- 
tributed to cosmic rays hitting the base of the detec- 
tor array were removed from the XQC spectrum, and 
these cuts also contribute to the dead time of the sys- 
tem. Specifically, events that arrived too close together 
for the calorimeters to accurately measure distinct ener- 
gies were discarded. This criterion removed 12% of the 
observed events and the resulting loss of sensitivity was 
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FIG. 4: Top panel: The XQC energy spectrum from 0-4 
keV in 5 eV bins. This spectrum does not have non-linearity 
corrections applied (see Ref. [5l|]). so the calibration lines at 
3312 eV and 3590 eV appear slightly below their actual en- 
ergies. The cluster of counts to the left of each calibration 
peak result from X-rays passing through the HgTe layer and 
being absorbed in the Si components where up to 12% of the 
energy may then be trapped in metastable states. Bottom 
panel: The XQC energy spectrum from - 2.5 keV in 5 eV 
bins. This spectrum, combined with the over-saturation rate 
of 0.6 events per second with energies greater than 4000 eV, 
was used in our analysis. 



included in the dead time of the calorimeters. When a 
cosmic ray penetrates the silicon base of the detector ar- 
ray, the resulting temperature increase is expected to reg- 
ister as multiple, nearly simultaneous, low-energy events 
on nearby calorimeters. To remove these events from the 
spectrum, we cut out events that were part of either a 
pair of events in adjacent detectors or a trio of events 
in any of the detectors that arrived within 3 ms of each 
other and had energies less than 2.5 keV. This procedure 
was expected to remove more than 97% of the events that 
resulted from cosmic rays hitting the base of the array. 
Nearly all of the events attributed to heating from cosmic 
rays had energies less than 300 eV, and a high fraction 
of the observed low-energy events were included in this 
cut. For example, seventeen of the observed twenty- four 
events with energies less than 100 eV were removed. The 
expected loss of sensitivity due to events being falsely 
attributed to cosmic rays was included in the calculated 
dead time of the calorimeters. Once all the dead time is 
accounted for, the 150 seconds of data collection is equiv- 
alent to 100.7 seconds of observation with all thirty-four 
calorimeters operational. 



The XQC calorimeters are capable of detecting energy 
deposits that exceed 20 eV, but full sensitivity is not 
reached until the energy surpasses 36 eV, and for approx- 
imately half of the detection time, the detector's lower 
threshold was set to 120 eV. The calorimeters cannot re- 
solve energies above 4 keV, and the 2.5-4 keV spectrum is 
dominated by the detector's interior calibration source: 
a ring of 2/iCi 41 Ca that generates Ka and K/3 lines at 
3312 eV and 3590 eV, respectively. We refer the reader 
to Ref. 51] for a complete discussion of the calibration of 
the detector. These limitations restrict the useful portion 
of the XQC spectrum to 0.03-2.5 keV. This spectrum is 
shown in Fig. [H along with the full spectrum from 0-4 
keV. The XQC field of view was centered on a region of 
the sky known to have an enhanced X-ray background 
in the 100-300 eV range, possibly due to hot gas in the 
halo, and this surge in counts can be seen in Fig. [4] In 
addition to the information present in this spectrum, we 
know that the XQC detector observed an average over- 
saturation event rate of 0.6 per second. This corresponds 
to a total of 60 events that deposited more than 4000 eV 
in a calorimeter. In Section llVBI we describe how we use 
the observed spectrum between 29 eV and 2500 eV and 
the integrated over-saturation rate to constrain the total 
cross section for elastic scattering between dark matter 
particles and nucleons. 



III. DETECTING DARK MATTER 
A. Incidence of dark matter particles 

The expected flux of dark matter particles into the de- 
tector depends on the density of the dark matter halo in 
the Solar System. Unfortunately, the local dark matter 
density is unknown and the range of theoretical predic- 
tions is wide. By constructing numerous models of our 
galaxy with various dark matter density profiles and halo 
characteristics, rejecting those models that contradict ob- 
servations, and finding the distribution of local dark mat- 
ter densities in the remaining viable models, Ref. [55| pre- 
dicted that the local dark matter density is between 0.3 
and 0.7 GeV cm -3 assuming that the dark matter halo 
is flattened, and the predicted local density decreases as 
the halo is taken to be more spherical. Another approach 
(56| used numerical simulations of galaxies similar to our 
own to find the dark matter density profile and then fit 
the profile parameters to Galactic observations, predict- 
ing a mean local dark matter density between 0.18 GeV 
cm~ 3 and 0.30 GeV cm' 3 . Given that it lies in the in- 
tersection of these two ranges, we use the standard value 
of 0.3 GeV cm~ 3 for the local dark matter density in 
our primary analysis. This assumption ignores the possi- 
ble presence of dark matter streams or minihalos, which 
do occur in numerical simulations [56| and could lead to 
local deviations from the mean dark matter density. 

We also assume that the velocities of the dark matter 
particles with respect to the halo are isotropic and have 
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a bounded Maxwellian distribution: the probability that 
a particle has a velocity within a differential volume in 
velocity-space centered around a given velocity v is 



P(v) 



iexp 




d v ifv< i>csc, 
if v > v csc . 



(2) 



where Vg is the dispersion velocity of the halo, v osc is the 
Galactic escape velocity at the Sun's position, and k is a 
normalization factor 1571: 



(™ 2 ) 3/2 



erf 



Vo 



vo 



(3) 



Numerical simulations indicate that dark matter particle 
velocities may not have an isotropic Maxwellian distri- 
bution [56| . Ref. [58j examines how assuming a more 
complicated velocity distribution would alter the flux of 
dark matter particles into an Earth-based detector. 

Given the flat rotation curve of the spiral disk at the 
Sun's radius and beyond and assuming a spherical halo, 
the local dispersion speed vq is the maximum rotational 
velocity of the Galaxy v c [59|]. Reported values for the 
rotational speed include 222±20 km s _1 [gjj], 228±19 km 
s" 1 [H, 184 ± 8 km s" 1 || and 230 ± 30 km s" 1 |3|. 
Recent measurements of the Galaxy's angular velocity 
have yielded values of f2 ga i = 28 ± 2 km s _1 kpc -1 [64| 
and 32.8 ± 2 km s^ 1 kpc" 1 65]. If the Sun is located 
8.0 kpc from the Galactic center, these angular velocities 
correspond to tangential velocities 224 ± 16 km s _1 and 
262 ± 16 km s _1 respectively. We adopt v c = 220 ± 30 
km s _1 as a centrally conservative value for the Galaxy's 
circular velocity at the Sun's location. 

The final parameter we need to obtain the dark mat- 
ter's velocity distribution is the escape velocity in the So- 
lar vicinity. The largest observed stellar velocity at the 
Sun's radius in the Milky Way is 475 km s _1 , which es- 
tablishes a lower bound for the local escape velocity [6(| . 
Ref. [13] used the radial motion of Carney-Latham stars 
to determine that the escape velocity is between 450 and 
650 km s _1 to 90% confidence, and Ref. [68[ obtained a 
90% confidence interval of 498 to 608 km s _1 from ob- 
servations of high- velo city stars. A kinematic derivation 
of the escape velocity [59| gives 



1 + ln 



Rq 



(4) 



where Rq is the distance from the Sun to the center of the 
Galaxy, and i? ga i is radius of the Galaxy. Observations of 
other galaxies suggest that our galaxy extends to about 
100 kpc and observations of Galactic satellites indi- 
cate that the Galaxy's flat rotation curve extends to at 
least 110 kpc (63[. The commonly accepted value for the 
Solar radius is Rq = 8.0 kpc 69]. Recent measurements 
include i? = 7.9 ± 0.3 kpc [zjand R a = 8.01 ± 0.44 kpc 
[7l| . and a compilation of measurements over the past 
decade [7l[ yields an average value of Rq = 7.80 ± 0.33 



kpc. To estimate the escape velocity, we use 100 kpc as 
a conservative estimate of the Galactic radius and the 
standard value i?o = 8.0 kpc. These parameters, com- 
bined with v c = 220 km s _1 , predict an escape velocity 
of 584 km s _1 , which falls near the middle of the ranges 
proposed in Refs. [67l. [68j. 

The isotropic Maxwellian velocity distribution given 
by Eq. © specifies the dark matter particles' motion 
relative to the halo. However, we are interested in their 
motion relative to the XQC detector: i/observcd = ?dm — 
^detector where the latter two velocities are measured with 
respect to the halo. The velocity of the detector with 
respect to the halo has three components: the velocity of 
the Sun relative to halo, the velocity of the Earth with 
respect to the Sun, and the velocity of the detector with 
respect to the Earth. 

When discussing these velocities, it is useful to de- 
fine a Galactic Cartesian coordinate system. In Galactic 
coordinates, the Sun is located at the origin, and the 
xy-plane is defined by the Galactic disk. The x-axis 
points toward the center of the Galaxy, and the y-axis 
points in the direction of the Sun's tangential velocity 
as it revolves around the Galactic center. The z-axis 
points toward the north Galactic pole and is antiparal- 
lel to the angular momentum of the rotating disk. The 
motion of the Sun through the halo has two compo- 
nents. First, there is the Sun's rotational velocity as 
it orbits the Galactic center: v c in the y direction. Sec- 
ond, there is the motion of the Sun relative to the spiral 
disk [z|: v Q = (10.00 ± 0.36, 5.25 ±0.62, 7.17 ±0.38) 
km s _1 in Galactic Cartesian coordinates. When the 
Earth's motion through the Solar System during its an- 
nual orbit of the Sun is expressed in Galactic coor- 
dinates [57]], the resulting velocity at the time of the 
XQC experiment (7.3 days after the vernal equinox) is 
^Earth = (29.14, 5.330, -3.597) km s _1 . 

The final consideration is the velocity of the detector 
relative to the Earth. The maximum velocity attained by 
the XQC rocket was less than 1.2 km s _1 . This velocity is 
insignificant compared to the motion of the Sun relative 
to the halo. Moreover, the XQC detector collected data 
while the rocket rose and while it fell, and the average 
velocity of the rocket was only 0.104 km s _1 . Therefore, 
we neglect the motion of the rocket in the calculation of 
the dark matter wind. Combining the motion of the Sun 
and the Earth then gives the total velocity of the XQC 
detector with respect to the halo during the experiment 
in Galactic Cartesian coordinates: i?dctector = (39.14 ± 
0.36, 230.5 ± 30, 3.573 ± 0.38) km s _1 . Subtracting the 
velocity vector of the detector relative to the halo from 
the velocity vector of the dark matter relative to the halo 
gives the dark matter's velocity relative to the detector in 
Galactic coordinates. However, we want the dark matter 
particles' velocities in the coordinate frame defined by 
the detector, where the z-axis is the field-of-view vector. 
The XQC field of view was centered on I = 90°, b = +60° 
in Galactic latitude and longitude [5l[, so the rotation 
from Galactic coordinates to detector coordinates may 
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be described as a clockwise 30° rotation of the z-axis 
around the x-axis, which is taken to be the same in both 
coordinate systems. 

B. Dark Matter Interactions 

Calorimetry measures the kinetic energy transferred 
from the dark matter to the absorbing material without 
regard for the specific mechanism of the scattering or 
any other interactions. Consequently, the dark matter 
detection rate for a calorimeter depends only on the mass 
of the dark matter particle and the total cross section for 
elastic scattering between the dark matter particle and an 
atomic nucleus of mass number A, which is proportional 
to the cross section for dark matter interactions with a 
single nucleon (arj n )- The calorimeter measures the recoil 
energy of the target nucleus (mass ror), 

T7 ( 1 2 \ 2m T m dm a , 

-Erec = 7i m dmV dm 7 ; ^(1 - COS# C m), (5) 

\2 J (m T + m dla y 

where TO d m and i»d m are dark matter particle's mass and 
velocity prior to the collision in rest frame of the target 
nucleus and 9cm is the scattering angle in the center-of- 
mass frame. 

If the momentum transferred to the nucleus, q 2 = 
2mT-E r ecj is small enough that the corresponding de 
Broglie wavelength is larger than the radius R of the 
nucleus (qR <C H), then the scattering is coherent. In co- 
herent scattering, the scattering amplitudes for each in- 
dividual component in the conglomerate body are added 
prior to the calculation of the cross section, so the total 
cross section is proportional to the square of the mass 
number of the target nucleus. Including kinematic fac- 
tors [45|, [73| , the cross section for coherent scattering off 
a nucleus is given by 



a coh (A) = A 2 



TO red (DM,Nuc) 



TO rod (DM,n) 



(6) 



where m rc d(DM, Nuc) is the reduced mass of the nucleus 
and the dark matter particle, m rc j(DM, n) is the reduced 
mass of a nucleon and the dark matter particle, and A is 
the mass number of the nucleus. Coherent scattering is 
isotropic in the center-of-mass frame of the collision. 

Dark matter particles may be massive and fast-moving 
enough that the scattering is not completely coherent 
when the target nucleus is large [zij. When the scat- 
tering is incoherent, the dark matter particle "sees" the 
internal structure of the nucleus, and the cross section for 
scattering is reduced by a "form factor," which is a func- 
tion of the momentum transferred to the nucleus during 
the collision (q) and the nuclear radius (R): 



da 

an 



O"coh 
47T 



F 2 (q,R). 



(7) 



Since q depends on the recoil energy, which in turn de- 
pends on the scattering angle, incoherent scattering is 
not isotropic. 



In this discussion of coherence, we have neglected the 
possible effects of the dark matter particle's internal 
structure by assuming that <Jr> n is independent of recoil 
energy. If the dark matter particle is not point-like then 
<TDn decreases as the recoil momentum increases due to a 
loss of coherence within the dark matter particle. Inco- 
herence within the dark matter particle has observational 
consequences [75[ , but these effects depend on the size of 
the dark matter particle. To avoid restricting ourselves 
to a particular dark matter model, we assume that the 
dark matter particle is small enough that nucleon scat- 
tering is always coherent; when we discuss incoherence, 
we are referring to the effects of the nucleus's internal 
structure. 

According to the Born approximation, the form fac- 
tor for nuclear scattering defined in Eq. Q is the 
Fourier transform of the nuclear ground-state mass den- 
sity J57L [76[ . The most common choice for the form factor 
IzIEl is F 2 (q,R) = exp[-(gi? rms ) 2 /(3fr 2 )], where i? rms 
is the root-mean-square radius of the nucleus. For a solid 
sphere, R 2 ms — (3/5)i? 2 , so this form factor is equivalent 
to the form factor used in Ref. [521 ] . This form factor is 
an accurate approximation of the Fourier transform of 
a solid sphere for (qR)/h <, 2, but it grossly underesti- 
mates the reduction in er for larger values of q [571 ] . The 
maximum speed of a dark matter particle with respect 
to the XQC detector is ~ 800 km s _1 (escape velocity 
+ detector velocity), and at that speed, the maximum 
possible value of qR/h for a collision with a Hg nucleus 
(A = 200) is nearly ten for a 100 GeV dark matter parti- 
cle, and the maximum possible value of qR/h increases as 
the mass of the dark matter particle increases. Clearly, 
this approximation is not appropriate for a large portion 
of the dark matter parameter space probed by the XQC 
experiment. 

Furthermore, a solid sphere is not a very realistic model 
of the nucleus. A more accurate model of the nuclear 
mass density is p(r) = J d 3 r' po(r')pi (r — r'), where po is 
constant inside a radius Rq = R 2 — 5s 2 and zero beyond 
that radius and p\ — exp[— r 2 /(2s 2 )], where s is a "skin 
thickness" for the nucleus . The resulting form factor 
is 



F(q,R) = 3 



sm(qR/h) - (qR/h) cos(qR/h) 



(qR/h) 3 



x exp 



-^(qs/h) 2 



(8) 



We follow Ref. [57|] in setting the parameters in Eq. ©: 
s = 0.9 fm and 

R 2 = [(l.23A 1/3 - 0.6) 2 + 0.6317T 2 - 5s 2 ] fm 2 , (9) 

where A is the mass number of the target nucleus. 

Despite its simple analytic form, the form factor given 
by Eq. ((SJ is computationally costly to evaluate repeat- 
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edly. We use an approximation: 



F 



exp 



9(0.81 

(qR/h) 



(f) {H s 
(*)" 



t exp 



9 fm 

0.9 fm 
R 



iff 



< 2, 



iff >2. 



(10) 

The low-g approximation combines the standard approx- 
imation for the solid sphere with the factor accounting 
for the skin depth of the nucleus. The high-g approxi- 
mation was derived from the asymptotic form of the first 
spherical Bessel function and normalized so that the to- 
tal cross section is as close as possible to the exact result. 
The error in the total cross section due to the use of the 
approximation is less than 1% for nearly all dark mat- 
ter masses; the sole exception is m,j m ~ 10 — 100 GeV, 
and even then the error is less than 5%. Unless other- 
wise noted, we use this approximation for the form factor 
throughout this analysis. We also assume that the dark 
matter particle does not interact with nuclei in any way 
other than elastic scattering. 



IV. ANALYSIS OF XQC CONSTRAINTS 

To obtain an accurate description of the XQC experi- 
ment's ability to detect strongly interacting dark matter 
particles, we turned to Monte Carlo simulations. The 
Monte Carlo code we wrote to analyze the XQC experi- 
ment simulates a dark matter particle's journey through 
the atmosphere to the XQC detector chamber, its path 
through the detector chamber to a calorimeter, and its in- 
teraction with the sensitive components of the calorime- 
ter. This latter portion of the code also records how much 
energy the particle deposits in the calorimeter through 
scattering. The results of several such simulations for the 
same set of dark matter properties may be used to pre- 
dict the likelihood that a given dark matter particle will 
deposit a particular amount of energy into the calorime- 
ter. These probabilities of various energy deposits predict 
the recoil-energy spectrum the XQC detector would ob- 
serve if the dark matter particles have a given mass and 
nucleon-scattering cross section. This simulated spec- 
trum may then be compared to the XQC data to find 
which dark matter parameters are excluded by the XQC 
experiment. 



A. Generating Simulated Energy-Recoil Spectra 

The basic subroutine in our Monte Carlo algorithm 
is the step procedure. The step procedure begins with 
a particle with a certain velocity vector and position in 
a given material and moves the particle a certain dis- 
tance in the material, returning its new position and ve- 
locity. The step procedure also determines whether or 
not a scattering event occurred during the particle's trek 



and updates the velocity accordingly. The number of ex- 
pected collisions in a step of length I through a material 
with target number density n is n x er tot x I, where <7t t 
is the total scattering cross section obtained by integrat- 
ing Eq. ([7]) over the scattering angle, or equivalently, the 
recoil momentum q: 



Ctot 



Ccoh 



9m a > 



F 2 (q,R) dq 2 



(11) 



where gwx is the maximum possible recoil momentum. 
The step length I is chosen so that it is at most a tenth 
of the mean free path through the material, so the num- 
ber of expected collisions is less than one and represents 
the probability of a collision. After each step, a ran- 
dom number between zero and one is generated using 
the "Mersenne Twister" (MT) algorithm [zi| and if that 
random number is less than the probability of a collision, 
the particle's energy and trajectory are updated. First, a 
recoil momentum is selected according to the probability 
distribution P(q 2 ) = F 2 (q, R)a coh /(q 2 na ^a to t), where the 
exact form factor is used for qR/h > 2 so that the oscilla- 
tory nature of the form factor is not lost. The recoil mo- 
mentum determines the recoil energy and the scattering 
angle in the center-of-mass frame through Eq. ([5]). The 
scattering is axisymmetric around the scattering axis, so 
the azimuthal angle is assigned a random value between 
and 2ir. The scattering angles are used to update the 
particle's trajectory, and its speed is decreased in accor- 
dance with the kinetic energy transferred to the target 
nucleus. The step subroutine repeats until the particle 
exits the simulation, or its kinetic energy falls below 0.1 
eV, or the energy deposited in the calorimeter exceeds 
the saturation point of 4000 eV. 

Our simulation treats the atmosphere as a 4.6x4.6 cm 
square column with periodic boundary conditions, the 
bottom face of which covers the top of the conical detec- 
tor chamber described in SectionQll This implementation 
assumes that for every particle that exits one side of the 
column, there is a particle that enters the column from 
the opposite side with the same velocity. The infinite 
extent of the atmosphere and its translational invariance 
makes this assumption reasonable. The atmosphere col- 
umn extends to an altitude 1000 km; increasing the at- 
mosphere height beyond 1000 km has a negligible effect 
on the total number of collisions in the atmosphere. The 
simulation begins with a dark matter particle at the top 
of the atmosphere column at a random initial position 
on the 4.6x4.6 cm square. Its initial velocity with re- 
spect to the dark matter halo is selected according to the 
isotropic Maxwellian velocity distribution function given 
by Eq. ([2]), and then the velocity relative to the detector 
is found via the procedure described in Section IIII Al 

The dark matter particle's path from the top of the at- 
mosphere to the detector is modeled using the step pro- 
cedure described above. The simulation of the particle's 
interaction with the atmosphere ends if the particle's al- 
titude exceeds 1000 km or if the particle falls below the 
height of the XQC rocket. We use the time-averaged alti- 
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tude (201.747 km) as the constant altitude of the rocket. 
We made this simplification because it allows us to ig- 
nore the periodic inactivity of each calorimeter and treat 
the detector as thirty-four calorimeters that are active for 
100.7 seconds. When the dark matter particle hits the 
rocket, its path through the five filter layers is also mod- 
eled using the step procedure, as is its path through the 
calorimeters. In addition to being smaller than the mean 
free path, the step length is chosen so that the particle's 
position relative to the boundaries of the detectors is ac- 
curately modeled. The simulation ends when the dark 
matter particle's random-walk trajectory takes it out of 
the detector chamber. As mentioned in Section |TT1 the 
calorimeter detects the sum of all the recoil energies if 
the dark matter particle is scattered multiple times. 

When the dark matter particle is unlikely to experience 
more than one collision in the calorimeter, this simula- 
tion is far more detailed than is required to accurately 
predict the energy deposited by the dark matter parti- 
cle. This is the case for the lightest (mdm < 10 2 GeV) 
and weakest-interacting (oDn < 10~ 26 cm 2 ) dark mat- 
ter particles that the XQC calorimeters are capable of 
detecting. Since the lightest dark matter particles are 
also the most numerous, many Monte Carlo trials are 
required to sample all the possible outcomes of a dark 
matter particle's encounter with the detector. The sim- 
ulation described above is too computationally intensive 
to run that many trials, so we used a faster and simpler 
simulation to model the interactions of these dark matter 
particles. This simulation assumes that the particle will 
experience at most one collision in the atmosphere and 
at most one collision in each filter layer and each layer of 
the calorimeter. The simulation ends if the probability 
of two scattering events in either the atmosphere or any 
of the filter layers exceeds 0.1. Instead of tracking the 
dark matter particle's path through the atmosphere, the 
total overburden for the atmosphere is used to determine 
the probability that the dark matter particle scatters in 
the atmosphere, and the particle only reaches the detec- 
tor if its velocity vector points toward the detector after 
the one allowed scattering event. Also, instead of the 
small step lengths required to accurately model the ran- 
dom walk of a strongly interacting particle, each layer is 
crossed with a single step. These simplifications reduce 
the runtime of the simulation by a factor of 100, making 
it possible to run 10 10 trials in less than one day. 



B. Comparing the Simulations to the XQC Data 

In order to compare the probability spectra produced 
by our Monte Carlo routine to the results of the XQC 
experiment, we must multiply the probabilities by the 
number of dark matter particles that are encountered by 
the initial surface of the Monte Carlo routine. When the 
initial velocity of the dark matter particle is chosen, the 
initial velocity may point toward or away from the detec- 
tor; in the latter case, the trial ends immediately. Con- 
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FIG. 5: Simulated event spectra for dark matter particles with 
masses of 1, 10 and 100 GeV and a total nucleon-scattering 
cross section of 10 -27 ' 3 cm 2 . In addition to the events de- 
picted in these spectra, the simulations predict 1300 ± 160 
events with energies greater than 4000 eV when mdm = 10 
GeV and 10, 000 ± 1200 such events when m dm = 100 GeV. 
The histogram represents the XQC observations. 



sequently, the Monte Carlo probability that the particle 
deposits no energy in the calorimeter already includes the 
probability that the dark matter particle does not have a 
halo trajectory that takes it into the atmosphere. There- 
fore, the probabilities resulting from the Monte Carlo 
routine should be multiplied by the number of particles 
in the volume swept out by the initial 4.6x4.6 cm 2 square 
surface during the 100.7/(i?) seconds of observation time, 
where / is the fraction of the observing time that the 
XQC detector was sensitive to deposits of energy E. For 
energies between 36 and 88 eV, / is 0.5083, and the value 
of / increases to one over energies between 88 and 128 eV. 
The detector was also slightly sensitive to lower energies: 
between 29 and 35 eV, / increases from 0.3815 to 0.5083. 
The normal of the initial surface points along the detec- 
tor's field of view, and the surface moves with the detec- 
tor; using the detector velocity given in Section UlI Al the 
number of dark matter particles encountered by the ini- 
tial surface is N Aai = / x (pdm/"Mm) x [(2.5 ±0.3) x 10 10 
cm 3 ], where pdm is the local dark matter density. 

The simulated event spectra produced by our Monte 
Carlo routine indicate that particles with masses less 
than 1 GeV very rarely deposit more than 100 eV inside 
the XQC calorimeters. Conversely, particles with masses 
greater than 100 GeV nearly always deposit more than 
4000 eV when they interact with the XQC calorimeters, 
so constraints on ctdii for these rndm values arise from the 



9 




Energy (eV) 

FIG. 6: Simulated event spectra for 10-GeV dark mat- 
ter particles with total nucleon-scattering cross sections of 
10~ 21 ' 6 , 10~ 27 ' 3 and 10~ 283 cm 2 . In addition to the events 
depicted in these spectra, the simulations predict 140 ± 37 
events with energies greater than 4000 eV when an n = 10 -21 ' 6 
cm 2 , 1300 ± 160 such events when ctdh = 10~ 27,3 cm 2 , and 
120±15 such events when arm = 10 -28 ' 3 cm 2 . The histogram 
represents the XQC observations. 



over-saturation (E > 4000 eV) event rate. Fig. [5] shows 
simulated spectra for three m^m values that lie between 
these two extremes, along with a histogram that depicts 
the XQC observations. Given an initial velocity of 300 
km s _1 relative to the XQC detector, a 1-GeV particle 
can only deposit up to 66 eV in a single collision with an 
Si nucleus, so the spectrum for these particles is confined 
to very low energies. Meanwhile, a 10-GeV particle and 
a 100-GeV particle with the same initial velocity can de- 
posit up to 900 eV and 44,000 eV, respectively, in a single 
collision with an Hg nucleus. In fact, ignoring any loss 
of coherence, all recoil energies between and 44,000 eV 
are equally likely during a collision between an Hg nu- 
cleus and a 100-GeV dark matter particle. That's why 
the mdm = 100 GeV spectrum in Fig.[5]is flat below 2500 
eV and why the simulations predict 10,000 events with 
energies greater than 4000 eV for this value of mdm and 

CTDn- 

Fig. [6] shows how changing the total cross section for 
elastic scattering off a nucleon affects the simulated spec- 
tra generated by our Monte Carlo routine for a single 
dark matter particle mass (rndm = 10 GeV). We see that 
increasing ctdh from 10~ 28 ' 3 cm 2 to 10~ 3 cm 2 increases 
all of the counts by a factor of ten but leaves the ba- 
sic shape of the spectrum unchanged. For much larger 
values of OD n , however, the particle loses a considerable 
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TABLE I: The binned XQC results used for comparison with 
our Monte Carlo simulations. 



amount of its energy while traveling through the atmo- 
sphere. Consequently, high-energy recoil events become 
less frequent, as shown by the spectrum for cr Dn = 10~ 216 
cm 2 . For larger values of ODn, too much energy is lost in 
the atmosphere for the particle to be detectable by the 
XQC experiment. 

When comparing the simulated measurements to the 
XQC data, we group the events into the thirteen energy 
bins given in Table HI We generally use large bins be- 
cause it reduces the fractional error in the probabilities 
generated by our Monte Carlo routine by increasing the 
probability of each bin: Spi jpi = 1/ \/pit, where t is the 
number of trials and pi is the probability of an energy 
deposit in the ith bin. Given that the number of trials is 
limited by runtime constraints, increasing the bin size is 
often the only way to obtain bin probabilities with bpi jpi 
values much less than one. When choosing our binning 
scheme, we attempted to maximize bin size while pre- 
serving as many features of the observed spectrum as 
possible. We also grouped all energies for which / ^ 1 
into two bins; we ignore the variation in / within these 
bins and set / = 0.3815 in the lowest-energy bin and 
/ = 0.5083 in the next-to- lowest bin. 

Unfortunately, we do not know the number of X-ray 
events in any of the bins listed in Table |U We considered 
using a model to subtract off the X-ray background but, 
given any model's questionable accuracy, we decided not 
to use it in our analysis. Our ignorance of the X-ray 
background forces us to treat the number of observed 
counts in each bin as an upper limit on the number of 
dark matter events in that energy range. Consequently, 
we define a parameter X 2 that measures the extent of 
the discrepancy between the simulated results for a given 
mdm and tiDn and the XQC observations while ignoring 
bins in which the observed event count exceeds the pre- 
dicted contribution from dark matter: 

*=# of Bins 2 
X 2 ^ Yl E with , (12) 

i=l ^ 1 ' 

where Ei — iVdm x pi is the number of counts in the ith 
bin predicted by the Monte Carlo simulation and Ui is 
the number of observed counts in the same bin. We use 
a second Monte Carlo routine to determine how likely it 
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is that a set of observations would give a value of X 2 as 
large or larger than the one derived from the XQC data 
given a mean signal described by the set of Ei derived 
from the simulation Monte Carlo. 

In the comparison Monte Carlo, a trial begins by gen- 
erating a new set of Ei by sampling the error distribu- 
tions of iVdm and pi. The distribution of iVdm values is 
assumed to be Gaussian with the mean and standard de- 
viation given above. The probability pi is derived from 
Pi x t events in the simulation Monte Carlo (recall that t 
is the number of trials), so a new value for pi is generated 
by sampling a Poisson distribution with a mean of pi x t 
and dividing the resulting number by t. Once a new set 
of Ei has been found, the routine generates a simulated 
number of observed counts for each bin according to a 
Poisson distribution with a mean of JSj . The value of the 
X 2 parameter for the new Ei and Ui is computed and 
compared to the value for the original Ei and the XQC 
observations, X 2 xqc- The number of trials needed to ac- 
curately measure the probability P(X) that X 2 > X 2 xqc 
is determined by requiring that the variation in the mean 
value of X 2 over ten Monte Carlo simulations does not 
exceed (100-C)%, where C% is the desired confidence 
level and that the range P(X)±(5x the variation in P(X)) 
does not contain (100 - C)/100. 



V. RESULTS AND DISCUSSION 

The XQC experiment rules out the enclosed region in 
(mdmjO'Dn) parameter space shown in Fig. [Jj The over- 
burden from the atmosphere and the filtering layers as- 
sures that there will be a limit to how strongly a dark 
matter particle can interact with baryons and still reach 
the XQC calorimeters; this overburden is responsible for 
the top edge of the exclusion region. Conversely, if aun is 
too small, the dark matter particles will pass through the 
calorimeters without interacting. The low-energy thresh- 
old of the XQC calorimeters places a lower bound on the 
excluded dark matter particle masses; if mdm is too small, 
then the recoil energies are undetectable. On the other 
side of the mass range, the XQC detector is not sensitive 
to mdm *J 10 5 GeV because the number density of such 
massive dark matter particles is too small for the XQC 
experiment to detect. 

The exclusion region shown in Fig.[7Jhas a complicated 
shape, but its features are readily explicable. As mdm in- 
creases, the range of excluded arm values shifts to lower 
values and then moves up again. The downward shift for 
mdm between 0.1 GeV and 100 GeV is due to the effects of 
coherent nuclear scattering. Since (T co h increases with in- 
creasing mdm for fixed ctdii, a 100-GeV particle interacts 
more strongly in the atmosphere and in the detector than 
a 1-GeV particle with the same <7Dn- Consequently, both 
the upper and lower boundaries of the excluded region 
decrease with increasing mass for md m & 100 GeV. The 
scattering of dark matter particles with larger masses is 
incoherent, and the form factor discussed in Section UlI Bl 
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FIG. 7: The region of dark matter parameter space excluded 
by the XQC experiment; o~D n is the total cross section for scat- 
tering off a nucleon and md m is the mass of the dark matter 
particle. This exclusion region follows from the assumption 
that the local dark matter density is 0.3 GeV cm -3 and that 
all of the dark matter shares the same value of or> n . 



causes er tot to decrease as mass increases for fixed (7Dn- 
Moreover, particles that are more massive than the target 
nuclei have straighter trajectories than lighter dark mat- 
ter particles due to smaller scattering angles in the de- 
tector rest frame. The loss of coherence also contributes 
because incoherent scattering makes small scattering an- 
gles more probable. A straight trajectory is shorter than 
a random walk, so the more massive particles interact 
less in the atmosphere and the detector than the more 
easily-deflected lighter particles. Due to both of these 
effects, the upper and lower boundaries of the exclusion 
region increase with increasing mdm for m,j m > 100 GeV. 

The lower left corner of the exclusion region also has 
two interesting features. First, the lower bound on the 
excluded value of ctdii decreases sharply as mdm increases 
from 0.1 GeV to 0.5 GeV. A dark matter particle with the 
maximum possible velocity with respect to the detector 
(800 km s^ 1 ) must have a mass greater than 0.24 GeV 
to be capable of depositing 29 eV in the calorimeter in 
a single collision. Lighter particles are only detectable 
if they scatter multiple times inside the calorimeter, and 
multiple scatters require a higher value of ctdii- Since 
their analysis does not allow multiple collisions, the XQC 
exclusion region found in Ref. [52| does not extend to 
masses lower than 0.3 GeV for any value of con- Second, 
there is a kink in the lower boundary at m,d m = 10 GeV; 
the constraint on tron is not as strong for this mass. The 
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FIG. 8: The region of dark matter parameter space excluded 
to 90% confidence by the XQC experiment for several values of 
the local density of dark matter with a total nucleon scattering 
cross section ODn and mass nid m . The four densities shown 
are 0.3 GeV cm -3 (solid line), 0.15 GeV cm -3 (long dashed 
line), 0.075 GeV cm" 3 (short dashed line) and 0.03 GeV cm" 3 
(dotted line). 



simulated spectra produced by our Monte Carlo routine 
for mdm = 10 GeV and ajy n 10" 25 cm 2 reveal that the 
particle is most likely to deposit between 100 and 600 
eV, as exemplified by the spectra depicted in Fig. [5] The 
background in this energy range is very high, so the XQC 
constraints are not as strict at these energies. 

Altering the local density of dark matter that strongly 
interacts with baryons changes the exclusion region. 
Fig. [8] shows the 90%-confidence exclusion regions for 
four values of the local density of dark matter particles 
with strong baryon interactions: 0.3 GeV cm" 3 (solid 
line), 0.15 GeV cm" 3 (long dashed line), 0.075 GeV cm" 3 
(short dashed line) and 0.03 GeV cm" 3 (dotted line). 
These different local densities could arise due to varia- 
tions in the local dark matter density due to mini-halos 
or streams. They also describe models where the dark 
matter does not consist of a single particle species and 
the dark matter that strongly interacts with baryons is 
a fraction fd of the local dark matter. In that case, the 
four exclusion regions in Fig. [H] correspond to f c i =1, 0.5, 
0.25, and 0.1. 

Fig. [H] indicates that the top and left boundaries of the 
XQC exclusion region are not highly sensitive to the dark 
matter density. In particular, the upper left corner of the 
exclusion region (0.01 < Wdm < 0.1 GeV) is nearly un- 
affected by lowering the dark matter density. This con- 



sistency indicates our Monte Carlo-generated exclusion 
region is smaller than the true exclusion region in this 
corner. If the dark matter is light (mdm ^ 0.1 GeV), 
then the number of dark matter particles encountered by 
the XQC detector is very large (A dm > 7 x 10 10 ). As 
previously mentioned, the upper left corner of the XQC 
exclusion region results from multiple scattering events, 
so the simpler version of our Monte Carlo code described 
in Section HV Al is not applicable. Consequently, it is not 
possible to run more than 10 9 trials in a week, so each 
scattering event in the simulation corresponds to more 
than one scattering event in the detector for all the den- 
sities shown in Fig. [8l Therefore, decreasing the density 
does not change the result. If it were possible to run 10 11 
trials, then the upper left corner of the exclusion region 
would expand and differences between the different den- 
sity contours would emerge. Since the upper left corner 
of the XQC exclusion region is already ruled out by astro- 
physical constraints (see Fig. [5]) , we have not invested in 
the computational time necessary to expand this corner. 

The upper boundary of the exclusion region is also not 
greatly affected by decreasing the particle density, even 
when TVdm is small enough that the Monte Carlo routine 
is capable of running more than Ad m trials (mdm > 100 
GeV). This robustness indicates that the overburden of 
the XQC experiment effectively prevents all dark matter 
particles with ctdii values greater than the upper bound 
of the exclusion region from reaching the detector, so 
it does not matter how many particles are encountered. 
Finally, as discussed previously, the lower portion of the 
exclusion region's left boundary (cr Dll < 10" 23 cm 2 ) is 
set by the energy threshold for detection and is therefore 
independent of Ad m . 

Examining the features of the excluded region allows 
us to predict how the region may be expanded by a fu- 
ture XQC- like experiment. Decreasing the overburden by 
either increasing the rocket's altitude or reducing the fil- 
tering will push the top boundary of the excluded region 
upwards. Decreasing the energy detection threshold will 
extend the excluded region to lower masses. It may also 
extend the exclusion region to higher values of crDn for all 
masses since strongly interacting particles lose much of 
their energy in the atmosphere and arrive at the calorime- 
ter with too little energy to produce a detectable sig- 
nal. Increasing the size or number of calorimeters would 
increase the sensitivity and extend the excluded region 
to lower values of a^ n . Finally, increasing the observa- 
tion time would increase Ad m , and that would extend the 
right and bottom boundaries of the excluded region. 



VI. CONCLUSION 

The X-ray Quantum Calorimetry (XQC) experiment is 
a powerful detector of dark matter that interacts strongly 
with baryons due to its high altitude and minimal shield- 
ing. The XQC measurements rule out a large range of 
hitherto unconstrained dark matter masses and scatter- 
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FIG. 9: Plot of the scattering cross section for dark matter particles and nucleons (udi) versus dark matter particle mass 
(mdm) showing the new XQC limits along with other current experimental limits. The red XQC exclusion region is the same 
as shown in Fig. [7] and the other experiments are discussed in the text. The dark gray region shows the maximal range of 
dark matter self-interaction cross section consistent with the strongly self-interacting dark matter model of structure formation 
[37l . l3Sj . The square marks the value of the scattering cross section for neutron-nucleon interactions. 



ing cross sections. The excluded range was first derived 
in Refs. [H, |48| based on rough analytic estimates. In 
this paper, we have improved upon these results using 
detailed Monte Carlo simulations to predict how a dark 
matter particle of a given mass and cross section for nu- 
cleoli scattering would interact with the XQC calorime- 
ters. Unlike Ref. [52] ]. our analysis includes the atmo- 
sphere and the shielding of the detector, so our result in- 
cludes the upper limit on excluded CTDn values, which had 
not yet been accurately determined. Our simulation also 
models the internal geometry of the XQC detector and 
the random walk of particles through it, which is not pos- 
sible using the analytical approaches of Refs. [38l. |48L152| . 

The resulting exclusion region is significantly different 
than its analytical predecessors. When multiple scatter- 
ings are included, the XQC experiment is sensitive to 
dark matter particles with masses below 0.3 GeV and 
cross sections for nucleon scattering between 10 -24 and 
KT 20 cm 2 . Unlike Ref. we find that the XQC ex- 
clusion region does not include a^ n < 10~ 29 cm 2 for 
dark matter masses less than 10 GeV. Ref. [52| obtained 
a more restrictive upper bound because they assumed a 
specific X-ray background while we treat all events as po- 
tential dark matter interactions. At higher masses, the 
lower boundar y o f our exclusion region is much higher 
than in Refs. [38l [48j because they over-estimated the 
XQC sensitivity by assuming coherent scattering. It also 



appears that Refs. [H, [48| underestimated the atmo- 
spheric and shielding overburden for the XQC detector 
because our exclusion region does not extend to values of 
CTDn as large as those included in their exclusion region. 
We also assume a lower local dark matter density than 
Refs. [HIH (°- 3 instead of 0.4 GeV cm" 3 ), so some of 
the shrinkage of the exclusion region may be attributed 
to the reduction in the assumed number density of dark 
matter particles. 

Fig.[5]shows how the XQC exclusion region depicted in 
Fig. [7] complements the exclusion regions from other ex- 
periments that are sensitive to similar values of CTo n and 
TOdm- For a summary of some of the other experimen- 
tal constraints as of 1994, see Ref. (46[. The constraints 
to £7 Dn from Pioneer 11 Q, Skylab (ph and IMP7/8 
[82| were interpreted by Refs. [H, I4H |48||. There have 
been two balloon-borne searches for dark matter, the 
IMAX experiment [H, H3| and the Rich, Rocchia & Spiro 
(RRS) [83| experiment. Although under grou nd detectors 
are designed to detect WIMPs, DAMA [H, H does ex- 
clude CTDn values within the range of interest, and relevant 
constraints may be derived from Edelweiss (EDEL) and 
CDMS [H,[I3- 

All of the exclusion regions shown in Fig. [9] were de- 
rived assuming that all the dark matter is strongly in- 
teracting. A local dark matter density of 0.4 GeV cm -3 
was assumed in the analysis of the exclusion regions from 
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Pioneer 11, Skylab and the RRS experiment, while all 
the other exclusion regions were derived assuming a lo- 
cal dark matter density of 0.3 GeV cm~ 3 . Furthermore, 
the derivations of all the shown exclusion regions other 
than the XQC region and the EDEL+CDMS region as- 
sume that the scattering between dark matter particles 
and nuclei is coherent. Therefore, these exclusion regions 
are likely too broad because they over-estimate the cross 
section for nuclear scattering. A comparison of the XQC 
exclusion region reported in Refs. [la, El] and our exclu- 
sion region indicates that assuming coherent scattering 
extends the exclusion region for mam > 1000 GeV to 
CTDn values that are roughly Ax smaller than the lower 
boundary of our exclusion region, where A is the mass 
number of the largest target nucleus. 

Fig. [5] also shows the bound on <jdii from the CMB 
and large-scale structure (LSS) obtained when one as- 
sumes prior knowledge of the Hubble constant Hq and 
the cosmic baryon fraction (from BBN) ■ This bound 
is nominally stronger than the bound from disk stabil- 
ity [45l ] , but it is less direct in that it requires combining 
different measurements and depends on the cosmologi- 
cal model; consequently, we show both bounds in Fig. [9] 
Measurements of primordial element abundances give an 
upper limit of cr Dn /m dm <J 4 x 10 -16 cm 2 GeV -1 (i^ . 
Since this upper bound lies well beyond the upper bound 
from disk stability, we do not include it in Fig. [5] We 
also do not display the constraints from cosmic rays (49| 
because they are derived from inelastic interactions that 
are model-dependent. 

As shown in Fig. [51 the XQC experiment rules out 
a wide region of (mdm, CDn) parameter space that was 
not probed by prior dark matter searches. Of partic- 
ular interest is the darkly shaded range of <7Dn values 
that corresponds to the maximal range of dark matter 
self-interaction cross sections consistent with the strongly 
self-interacting dark matter model of structure formation 
[37l . [Hj]. If the dark matter consists of exotic hadrons 
whose interactions with nucleons are comparable to their 
self-interactions, then aon for these particles would lie in 
or near the darkly shaded region in Fig. [5J Previous esti- 
mates of the XQC exclusion region [38, 48] indicated that 
the XQC experiment rules out all the darkly shaded aun 
values for 1 <, m dn i ^ 10 4 GeV. Our analysis reveals that 
this is not the case; portions of the darkly shaded region 
for mdm <2 20 GeV are not excluded by the XQC ex- 
periment, although they are ruled out by observations of 
LSS and the CMB. The mass-er combination correspond- 
ing to nucleon-neutron scattering (the square in Fig. [5]) 
lies within the exclusion region of the XQC experiment, 
and the only portion of the darkly shaded region that is 



unconstrained corresponds to dark matter masses smaller 
than 0.25 GeV. 

It is important to note, however, that the cross section 
for dark matter self-interactions need not be comparable 
to the cross section for nucleon scattering; a^n could dif- 
fer by a few orders of magnitude from the self-interaction 
cross section (as is the case for Q-balls). Furthermore, no 
interactions with baryons are required for self-interacting 
dark matter to resolve the tension between the collision- 
less dark matter model and observations of small-scale 
structure. 

Another XQC detector is scheduled to launch in the 
upcoming year. This experiment will have twice the ob- 
serving time of the XQC experiment used in this analy- 
sis. As discussed in Section [Vj increasing the observing 
time will extend the exclusion region to higher masses 
and weaker interactions. The future XQC experiment 
will also have a lower energy threshold (15 eV) and will 
maintain sensitivity to all energies above this threshold 
throughout the run. The increased sensitivity to low en- 
ergies will shift the lower (ctdii < 10~ 23 cm 2 ) left bound- 
ary of the exclusion region to lower masses. A lower en- 
ergy threshold of 15 eV will make the experiment sensi- 
tive to single recoil events involving dark matter particles 
more massive than 0.17 GeV, as discussed in Section [Vl 
Clearly, the next-generation XQC experiment will be an 
even more powerful probe of interactions between dark 
matter particles and baryons than its predecessor. 
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